Universal quantum criticality at the Mott- Anderson transition 
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We present a large N solution of a microscopic model describing the Mott-Anderson transition 
on a finite-coordination Bethe lattice. Our results demonstrate that strong spatial fluctuations, due 
to Anderson localization effects, dramatically modify the quantum critical behavior near disordered 
Mott transitions. The leading critical behavior of quasiparticle wavefunctions is shown to assume 
a universal form in the full range from weak to strong disorder, in contrast to disorder-driven non- 
Fermi liquid ("electronic Griffiths phase") behavior, which is found only in the strongly correlated 
regime. 
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Both Anderson (disorder-driven) and Mott 
(interaction-driven) routes to localization play sig- 
nificant roles in materials close to the metal-insulator 
transition (MIT) [T), but the full understanding of their 
interplay remains elusive. These effects are believed to 
be the cause of many puzzling features in systems rang- 
ing from doped semiconductors [2] and two-dimensional 
electron gases [3] , to broad families of complex oxides [4] . 

Early theories of the MIT [SJ focused on the stability 
of conventional metals to introducing weak disorder, but 
these Fermi-liquid approaches proved unable to describe 
strong correlation phenomena which took center stage 
in recent years. Indeed, the last two decades have seen 
significant advances in our understanding of "Mottness" , 
not least because of the development of dynamical mean- 
field theory (DMFT) [6], which has been very successful 
in quantitatively explaining many aspects of strong cor- 
relation. 

Because in its simplest form DMFT does not capture 
Anderson localization effects, several recent refinements 
were introduced, which proved capable of capturing both 
the Mott and the Anderson mechanisms. The conceptu- 
ally simplest [7] such approach - the typical-medium the- 
ory (TMT) - provided the first self-consistent description 
of the Mott-Anderson transition, and offered some insight 
into its critical regime. For weak to moderate disorder, 
this theory found a transition closely resembling the clean 
Mott point, while only at stronger disorder Anderson lo- 
calization modified the critical behavior. However, close 
scrutiny [5] revealed that some of these findings may be 
artifacts of the neglect of spatial fluctuations in this for- 
mulation. 

An alternative but technically more challenging ap- 
proach to Anderson-Mott localization was dubbed "Sta- 
tistical DMFT" (staiDMFT) [5]. Here, only the strong 
correlations are treated in a self-consistent DMFT fash- 
ion, while disorder fluctuations are treated by a (nu- 
merically) exact computational scheme. While certainly 
much more reliable than TMT-DMFT, so far this method 



was utilized only in a handful of theoretical studies of the 
Mott-Anderson transition [10] [11] , and the precise form 
of quantum criticality has never been explored in detail. 

In this letter, we present the first precise and de- 
tailed study of the quantum critical behavior of the Mott- 
Anderson transition in a half-filled Bethe lattice, within 
the framework of statDMFT. We address the following 
physical questions, which we answered in a clear and reli- 
able fashion: (1) Are there two distinct types of quantum 
criticality in this model, as TMT-DMFT suggested, or do 
the fluctuation effects restore universality within the crit- 
ical regime? (2) How general is the disorder-driven non- 
Fermi liquid behavior (electronic Griffiths phase), and 
how does it relate to the relative strength of correlations 
and disorder? These results are obtained for a specific 
microscopic model, where our formulation proves exact 
in an appropriately defined large N limit. 

Charge transfer model and statDMFT - The charge 
transfer (CT) model has been used in the descrip- 
tion of many multi-band systems near the Mott transi- 
tion, including various oxides [T2], two-dimensional elec- 
tron gases near Wigner-Mott crystallization [T3], as well 
as doped semiconductors [14]. It consists of a two 
band model, where one band represents weakly corre- 
lated conduction electrons and the other one is a nar- 
row, strongly correlated band containing nearly Mott- 
localized electrons. In the clean case, the Mott-insulating 
phase can be approached at half-filling by increasing the 
charge-transfer gap, which reduces the hybridization be- 
tween the two bands. When the effective (correlation- 
renormalized) hybridization vanishes, the system under- 
goes a Mott metal-insulator transition. 

While many studies of the Mott-Anderson transitions 
concentrate on disordered single-band Hubbard models, 
we chose to focus on the CT model for two specific 
reasons. First, most real Mott systems are better de- 
scribed by the CT model than by the single-band Hub- 
bard model; the latter features unrealistic particle-hole 
symmetry which could modify [5] the critical behavior. 
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FIG. 1: (Color online) Phase diagram of the disordered CT 
model. E ct = —Ef is the CT energy, which plays the role of 
Hubbard's U. The dashed lines are guides to the eyes. 



Second, the simplest version of the CT model imposes the 
U = oo constraint (no double-occupancy) for the corre- 
lated band, a feature that considerably simplifies some 
technical aspects of our calculation, allowing a simple 
Gutzwiller-type variational solution of the model. This 
becomes formally exact in an appropriate large N limit, 
which is conveniently formulated within the standard 
slave-boson approach [15] , 

The Hamiltonian of our disordered CT model corre- 
sponds to that of the disordered Anderson lattice model, 
which is given by 

H = J2 K e < - A 5 ** - *] C U. + (Ef - m) ]T flfir 
+ v /2 ( c lfi« + ft,*,) +Uj2nmn fih (1) 



where c\ a (ci a ) creates (destroys) a conduction electron 
with spin a on site i, fj a and /j CT are the correspond- 
ing creation and annihilation operators for a localized 
/-electron with spin a on site i, Ufi a = j\ a ha is t ne num- 
ber operator for /-electrons, t is the hopping amplitude, 
Ef is the /-electron energy, U is the on-site repulsion 
between /-electrons, V is the hybridization between con- 
duction and /-electrons, and /i is the chemical potential. 
Disorder is introduced through the on-site energies for 
conduction electrons, which follow a distribution P(e), 
assumed to be Gaussian with zero mean and variance 
W 2 . Throughout this paper we use the half-bandwidth 
for conduction electrons as the unit of energy; the hy- 
bridization potential is chosen to be V = 0.5. 

Within the CT model, the average number of electrons 
per unit cell is equal to 1, which can be enforced by ad- 
justing the chemical potential and is written as 



1. 



(2) 



where nji = n + n gives the number of /-electrons 
on site i, n C i — n c ^ + n C i^ is the corresponding number 
operator for conduction electrons, with n C i a = c\ a Ci a , 
and the averages are taken over the distribution P(e). 

In this work we solve the above Hamiltonian using 
statDMFT [9 , in which the disordered lattice model is 
reduced to a self-consistent solution of an ensemble of 
single-impurity models located in the different sites j of 
the lattice. The corresponding local effective actions take 
the form: 
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+ A/i(r - r')] / JCT (r') + U / dTn m (T)n m (r), 
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where the superscript (i) indicates that site i has been 
removed from the lattice. The bath function, A/j(r— r'), 
is given by 



v 2 



iuj + fi — Ej — A c j (iuS) 



(4) 



where A C j(iuj) — t 2 X)fc=i G^(iu), k ^ i in the summa- 
tion, and z is the coordination number (chosen to be 3 in 
this work). G^(uj), the Green's function for conduction 
electrons, satisfies 



iu> + ii - e k - A ck (iuj) - $fe(iw) ' 



(5) 



where 



V 



■fi- Ef - T,f k {iu>) 



(6) 



and £/fc(«j) is the single-impurity self-energy, which is a 
solution of an action similar to that given in eq. (J3j>, but 
for site k (instead of jf). 

To solve the single- impurity problems of eq. ^ , we use 
the slave-boson (SB) technique in the U —> oo limit |15j . 
In this case, the impurity Green's function can be written 
as 



Gf k (iuj) 



ius - e fk - Z k Af k (iu) ' 



(7) 



where Z k is the local quasi-particle (QP) weight and Ef k 
is the renormalized /-electron energy. For technical de- 
tails on solving the statDMFT equations for a closely 
related model see Ref. [TfJ] , 

Order parameter and phase diagram - To determine 
the phase diagram within statDMFT, we examine the 
behavior of the typical value of the local density-of-states 
(LDOS) for conduction electrons at the Fermi energy 
(u> = 0). This order parameter [5] measures the degree 
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FIG. 2: (Color online) Typical (open symbols) and inverse 
of average (full symbols) values of the LDOS for conduction 
electrons close to the Fermi energy. Results are shown as a 
function of the disorder strength, W, for different values of 
the CT energy, E ct . The inset gives a zoom on the results 
close to the MIT. 

of localization of quasiparticle wavefunctions [7J [8] ■ It is 
given by 

Ptyp = exp{ (In pi (u = 0))}, (8) 

where (...) denotes the (arithmetic) average over disor- 
der and the most probable value is represented by the 
geometric average. Its vanishing directly indicates that 
wavefunctions are Anderson-localized. 

Fig. [I] presents the phase diagram we obtained for the 
disordered CT model. Here, we show the phase boundary 
between the metallic and the Mott insulating phases, as 
well as that between the metal and the Anderson insula- 
tor; the region where a metallic Griffiths phase is seen is 
also shown, close to the Mott insulator. For our model, 
the charge-transfer gap, E ct — —Et, controls the effec- 
tive strength of electron-electron interactions, i.e. plays 
the role of the Hubbard U. 

In dramatic contrast to TMT-DMFT jS] and also 
statDMFT [IT] results for disordered Hubbard models, 
here (see Fig. [TJ we find a surprisingly broad intermediate 
metallic regime, spanning the region where correlations 
and disorder are comparable. This result may indicate 
that strong correlation effects have, for charge-transfer 
models, a more pronounced tendency to suppress Ander- 
son localization than what was expected from Hubbard 
model studies. This finding is significant, because the CT 
model should be considered a more accurate representa- 
tion of many real materials than the simple single-band 
Hubbard model. 

Critical behavior near the disorder- driven transition - 
To obtain the phase boundary between the metallic and 
the Anderson insulating phases of Fig. [T] we examined 
the behavior of pt yp as a function of W, for different 
values of the CT energy; typical results are shown in 



FIG. 3: (Color online) Typical (open symbols) and inverse 
of average (full symbols) values of the LDOS for conduction 
electrons close to the Fermi energy. Results are shown as 
a function of the CT energy, E c t, for different values of the 
disorder strength W. 

Fig. [2] For all values of E ct considered, p typ is found 
to decrease with disorder, reflecting Anderson localiza- 
tion effects. As the interaction parameter E ct increases, 
we observe that the MIT shifts to larger values of dis- 
order. This behavior, which is caused by the tendency 
for disorder screening by correlation effects, is consistent 
with previous results obtained for the disordered Hub- 
bard model within TMT-DMFT 8 . In addition, we find 
that the inverse of p av = (pi(u = 0)) (see Fig. [2]) vanishes 
at exactly the same W value for which p typ vanishes, sim- 
ilarly as in statDMFT studies of disorder-driven MIT in 
doped Mott insulators [S]. 

Critical behavior near the Mott-like transition - The 
behavior of p typ as a function of E ct near interaction- 
driven (Mott-like) transitions is shown in Fig. [3] We 
find that, as disorder increases, the transition (given by 
Ptyp = 0) shifts to larger E ct , consistent with previous 
results for disordered Hubbard models [17] , reflecting 
the tendency of disorder to broaden the Hubbard bands. 

The behavior of our localization order parameter p typ 
shows very interesting behavior as the transition is ap- 
proached. It initially increases with the interaction pa- 
rameter E c t, reflecting the correlation-induced screening 
of disorder within the metallic phase, an effect already 
found in previous studies [T7J . Very close to the tran- 
sition, however, this tendency is reversed, and our order 
parameter p typ starts to decrease and eventually van- 
ishes, indicating localization of quasiparticle wavefunc- 
tions. Most remarkably, the critical behavior of both p typ 
and the inverse of p av — (pi(uj = 0)) vanish linearly at 
criticality - thus displaying precisely the same critical be- 
havior as for the disorder-driven transition. This result is 
in dramatic contrast to the prediction of simpler DMFT 
or TMT-DMFT treatments [8"1HT|. where p typ displays a 
finite jump at the disordered Mott transition. 
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The electronic Griffiths phase - The second order pa- 
rameter identified by statDMFT is the local quasiparticle 
(QP) weight Z [6j [9] ; its vanishing indicates the destruc- 
tion of the Fermi liquid by local moment formation (i.e. 
local Mott localization). In presence of disorder, this 
quantity displays spatial fluctuations, and must be char- 
acterized by an appropriate distribution function, P(Z), 
which can assume a power-law form P(Z) ~ Z Q_1 [9],I18|. 
When the corresponding exponent a becomes smaller 
than 1, this indicates [9] disorder-driven non- Fermi liq- 
uid behavior Q], dubbed the "electronic Griffiths phase" 
(EGP). 

We have carefully studied the behavior of P(Z) across 
the phase diagram, and we find that EGP behavior is 
not found close to disorder-driven transitions, but is re- 
stricted to the vicinity of the interaction-driven (Mott- 
like) transition (see Fig. 1). Similar behavior has also 
been observed for the two-dimensional Hubbard model 
within statDMFT [18]. We also found that a finite frac- 
tion of sites display Z — within the Mott-like insulator, 
confirming the formation of localized magnetic moments 
("Mott droplets") [18] . In contrast, all Z-s remain finite 
within the Anderson-like insulator, indicating the lack of 
local moment formation in this phase. 

Comparison with DMFT results - For further insight in 
the role of Anderson localization in our model, we com- 
pare our statDMFT results with those we obtained within 
conventional DMFT [T7], which can be easily adapted 
for our model. Here, the average Green's function for 
conduction electrons satisfies the following self-consistent 
equation 



G c {iu>) 



1 



■ + p-e k -t 2 G(iu)-<S> k {iuj) 



(9) 



Since conventional DMFT is unable to capture Ander- 
son localization effects but is able to describe Mott lo- 
calization, the comparison is only meaningful and inter- 
esting close to the interaction-driven transition. Within 
DMFT the typical and the average values of LDOS coin- 
cide, and in Fig.|4]we compare the DMFT results with the 
corresponding quantities obtained from our statDMFT 
calculation. Further from the transition, for small E ct , 
results obtained from both theories essentially coincide, 
displaying correlation-induced screening. Closer to the 
transition, p typ obtained from statDMFT displays a max- 
imum and starts to deviate from the DMFT result, and 
eventually vanishes much before DMFT would predict 
Mott localization (see Fig. 4). This clearly indicates 
that Anderson localization effects "kick in" only suffi- 
ciently close to Mott localization, and ultimately deter- 
mine the resulting quantum critical behavior. We stress 
here that its precise (linear) form, though, is qualitatively 
different [5j than in the non-interacting model, where 
Ptyp vanishes exponentially for the Bethe lattice model 
we consider. 
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FIG. 4: (Color online) StatDMFT typical (pt yP ) and inverse 
of average (l/p a v) and DMFT average (pdmft) values of the 
LDOS for conduction electrons close to the Fermi energy. 
Results are shown as a function of the CT energy, E ct , for 
(a) W = 0.25 and (b) W = 1.0. Note that both averages 
predicted by statDMFT start to deviate from DMFT results 
only within the narrow critical region close to the Mott-like 
transition. 



We reiterate that these disorder-induced localization 
effects arise only within a narrow critical region close to 
the disordered Mott insulator. Its size can be estimated 
by examining how the location of pt yp maximum evolves 
with disorder W . Remarkably, we find that the bound- 
ary of this critical region closely tracks the boundary of 
the EGP. This result makes perfect sense, as both An- 
derson localization effects and the EGP behavior reflect 
profound and dramatic spatial fluctuations of the rele- 
vant physical quantities. 

Conclusions - In summary, we carried out a detailed 
statDMFT study of the T = critical behavior at 
the Mott-Anderson transition for a disordered charge- 
transfer model. Remarkably, identical critical behavior 
of the localization order parameter p typ is found in the 
close vicinity of both interaction-driven (Mott-like) and 
disorder driven (Anderson-like) transitions. This result 
demonstrates that (spatial) fluctuation effects captured 
by our statDMFT, but not by simpler treatments, restore 
a degree of universality at this quantum critical point. In 
contrast, the thermodynamic behavior characterized by 
non-Fermi liquid features (electronic Griffiths phase) is 
found only near the Mott-like transition. Finite temper- 
ature properties as well as the role of inter-site correla- 
tions remain to be examined in more detailed, but these 
interesting questions remain a challenge for future work. 
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